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Abstract 

While the properties of the 125 GeV Higgs boson-like particle observed by the ATLAS and 
CMS collaborations are largely compatible with those predicted for the Standard Model state, 
significant deviations are present in some cases. We, therefore, test the viability of a Beyond 
the Standard Model scenario based on Supersymmetry, the CP- Violating Next-to-Minimal Su- 
persymmetric Standard Model, against the corresponding experimental observations. Namely, 
we identify possible model configurations in which one of its Higgs bosons is consistent with 
the LHC observation and evaluate the role of the explicit complex phases on both the mass 
and di-photon decay of such a Higgs boson. Through a detailed analysis of some benchmark 
points corresponding to each of these configurations, we highlight the impact of the CP-violating 
phases on the model predictions compared to the CP-conserving case. 
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1 Introduction 



In July 2012, the CMS and ATLAS experimental collaborations at the Large Hadron Collider 
(LHC) announced the observation of a new boson [H[2], consistent with a Higgs particle, the last 
undiscovered object in the Standard Model (SM). The initial results were based on data corre- 
sponding to integrated luminosities of 5.1 fb _1 taken at y/s = 7TeV and 5.3 fb -1 at 8TeV and the 
search was performed in six decay modes: H — > 77, ZZ, Z*y, WW, t + t~ and bb. A ~5cr excess of 
events with respect to the background was clearly observed in the first and second of these decay 
modes, while the remaining ones yielded exclusion limits well above the SM expectation. Both 
collaborations have since been regularly updating their findings [3l HI [5j O [7] , improving the mass 
and (so-called) 'signal strength' measurements. 

In these searches, the magnitude of a possible signal is characterised by the production cross 
section times the relevant Branching Ratios (BRs) relative to the SM expectations in a given channel 
X, denoted by R(X) = a/asM x BR(Y)/BRsm(Y) (i-e., the signal strength). According to the 
latest results released by the two collaborations after the collection of ~20 fb _1 of data El [7J, 
a broad resonance compatible with a 125 GeV signal is now also visible in the WW — >• 2l2v decay 
channel. The mass of the observed particle is still centered around 125 GeV but the measured 
values of its signal strength in different channels have changed considerably compared to the earlier 
results. These values now read 

R(H -> 77) 
R(H -> ZZ) 
R(H -> WW) 

at CMS, and 

R(H -> 77) 
R(H -»■ ZZ) 
R(H -> WW) 

at ATLAS. The bulk of the event rates comes from the gluon-gluon fusion channel [8]. Furthermore, 
the signal has also been corroborated by Tevatron analyses [9], covering the bb decay mode only, 
with the Higgs boson stemming from associated production with a W boson [8j. However, there 
the comparisons against the SM Higgs boson rates are biased by much larger experimental errors. 

If the current properties of the observed particle are confirmed after an analysis of the full 7 
and 8TeV data samples from the LHC, they will not only be a clear signature of a Higgs boson, 
but also a significant hint for possible physics beyond the SM. In fact, quite apart from noting that 
the current data are not entirely compatible with SM Higgs boson production rates, while the most 
significant LHC measurements point to a mass for the new resonance around 125 GeV the Tevatron 
excess in the bb channel points to a range between 115 GeV and 135 GeV. While the possibility that 
the SM Higgs boson state has any of such masses would be merely a coincidence (as its mass is a 
free parameter), in generic Supersymmetry (SUSY) models the mass of the lightest Higgs boson 
with SM-like behaviour is naturally confined to be less than 180 GeV or so [10]. The reason is that 
SUSY, in essence, relates trilinear Higgs boson and gauge couplings, so that the former are of the 
same size as the latter, in turn implying such a small Higgs boson mass value. Therefore, the new 
LHC results could well be perceived as being in favour of some low energy SUSY realisation. 

Several representations of the latter have recently been studied in connection with the aforemen- 
tioned LHC and Tevatron data, including the Minimal Super symmetric Standard Model (MSSM) 
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[11] (also the constrained version [12] of it, in fact), the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) p31 M, CES] , the E 6 -inspired Supersymmetric Standard Model (E 6 SSM) pi] and 
the (B-L) Supersymmetric Standard Model ((B-L)SSM) [17J. All of these scenarios can yield a 
SM-like Higgs boson with mass around 125 GeV and most of them can additionally explain the 
excesses in the signal strength measurements in the di-photon channel. 

Another approach to adopt in order to test the viability of SUSY solutions to the LHC Higgs 
boson data is to consider the possibility of having CP-violating (CPV) phases (for a general review 
of CP Violation, see [H]) in (some of) the SUSY parameters. These phases can substantially modify 
Higgs boson phenomenology in both the mass spectrum and production/decay rates at the LHC 
[191 l20l [2~T1 [22] , while at the same time providing a solution to electroweak baryogenesis |23] . In the 
context of the LHC, the impact of CPV phases was emphasised long ago in [241 [25] and revisited 
recently in [26] following the Higgs boson discovery. In all such papers though, CPV effects were 
studied in the case of the MSSM. 

In this paper, we consider the case of similar CPV effects in the NMSSM. In particular, we 
study the possibility to have Higgs boson signals with mass around 125 GeV in the CPV NMSSM 
which are in agreement with the aforementioned LHC data as well as the direct search constraints 
on sparticle masses from LEP and LHC. We also investigate the dependence of the feasible CPV 
NMSSM signals on the mass of the Higgs boson as well as its couplings to both the relevant particle 
and sparticle states entering the model spectrum, chiefly, through the decay of the former into a 77 
pair. We thus aim at a general understanding of how such observables are affected by the possible 
complex phases explicitly entering the Higgs sector of the next-to-minimal SUSY Lagrangian. 

The paper is organised as follows. In the next section we will briefly review the possible explicit 
CPV phases in the Higgs sector of the NMSSM. In Sect. [3] we will outline the independent CPV 
NMSSM parameters and the methodology adopted to confine our attention to the subset of them 
that can impinge on the LHC Higgs boson data. In the same section, we further investigate the 
possible numerical values of the complex parameters after performing scans of the low energy CPV 
NMSSM observables compatible with the LEP and LHC constraints on Higgs boson and SUSY 
masses. In Sect. [4] we present our results on the Higgs boson mass spectrum as well as signal rates 
in connection with the LHC. Finally, we conclude in Sect.[5j 

2 CPV phases in the Higgs sector of the NMSSM 

The CPV phases appearing in the Higgs potential of the NMSSM at tree-level [27] can be divided 
into three categories: 

1. 6 and ip: the spontaneous CPV phases of the vacuum expectation values (vevs) of the up- type 
Higgs doublet H u and the Higgs singlet S, respectively, with respect to the down-type Higgs 
doublet Ha; 

2. <p\ and 4> K : the phases of the Higgs boson trilinear couplings A and k; 

3. (pA x an d 4>A K '- the phases of the trilinear soft terms A\ and A K . 

As explained in [28[ [29] the phases in 3. above are determined by the minimisation conditions 
of the Higgs potential with respect to the three Higgs fields. Furthermore, assuming vanishing 
spontaneous CPV phases in 1., the only actual physical phases appearing in the tree- level Higgs 
potential are those in 2. as the difference (ft\ — <f) K . Beyond the Born approximation, the phases of the 
trilinear couplings At, Af, and A T also enter the Higgs sector through radiative corrections from the 
third generation squarks and stau (assuming negligible corrections from the first two generations). 
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Also, in the one-loop effective potential, each of the phases in 2. above contributes independently 
from the other. The complete one-loop Higgs mass matrix can be found in [28} l29j 130] . Here we 
only reproduce the tree-level Higgs as well as the sfermion mass matrices in Appendix A. 

The 5x5 Higgs mass matrix M. 2 N (after rotating out the Goldstone modes) is diagonalised with 
a unitary matrix O to yield five mass eigenstates as 

(<tP d ,<t>l,4>%,a,as) T = O {H u H 2 ,H z ,H^H b ) T , (1) 

where T Ai 2 N = diag(m|^ 1 ,m| f2 ,m^ , m^ 4 ,m|j g ) in order of increasing mass. For a non-zero 
value of any of the phases listed above, these mass eigenstates become CP-indefinite due to scalar- 
pseudoscalar mixing. Moreover, these CPV phases not only affect the masses of the Higgs states 
but also their decay rates (both widths and BRs), since the Higgs boson couplings to various 
particles are proportional to the elements of the unitary matrix O (see, e.g., [311 132] ). Additionally, 
modifications in the masses of light neutralinos, in particular, due to the CPV phases can also have 
an indirect yet non- negligible impact on the BRs of the Higgs boson into SM particles. 

The decay widths and BRs of the Higgs boson in the NMSSM with CPV phases can be calculated 
using the methodology implemented in [31]. Explicit expressions for Higgs boson couplings and 
widths in the CPV NMSSM can be found in [33], which follows the notation of [34]. These widths 
and BRs can then be used to obtain the signal strength of the 77 channel (also called reduced 
di-photon cross section), Rff t , defined as 

Rll = a(gg -» Hj) ^ BR(^ -» 77) 
H < a(gg^h SM ) BR(h SM -)■ 77) ' 

for a given Higgs boson, Hi. In terms of the reduced couplings Ci(X) (couplings of Hi with respect 
to those of the SM Higgs boson with the same mass), eq. ([2]) can be approximated by 

*g = [Qto)] 2 [Q(77)] 2 E B ^X^ X) ' (3) 
where X runs over all possible decay modes of the Higgs boson in the SM. 



3 Model parameters and methodology 

In light of the recent LHC discovery of a SM Higgs boson-like boson we scan the parameter space 
of the CPV Higgs sector of the NMSSM using a newly developed FORTRAN code. In our scans the 
LEP constraints on the model Higgs bosons are imposed in a modified fashion, i.e., they have to 
be satisfied by the scalar and pseudoscalar components of all the CP-mixed Higgs bosons. Also 
imposed are the constraints from the direct searches of the third generation squarks, stau and the 
light chargino at LEP. We should also point out that, in the CPC limit, the Higgs boson mass and 
BRs have been compared with those given by NMSSMTools [35] and have been found to differ from 
the latter by ~1% and ~5% at the most, respectively. Although no limits from 6-physics or from 
relic density measurements have been imposed we, nevertheless, confine ourselves to the regions of 
the parameter space which have been found to best comply with such constraints (see, e.g., [15]). 

We study the effects of the CPV phases described in the previous section on the mass and 
di-photon signal rate of a Higgs boson predicted by the model that is compatible with the Higgs 
boson discovery data from the LHC. In particular, we consider three most likely scenarios specific 
to the CPV NMSSM that comply with the latter. In our analyses, we assume minimal Supergravity 
(mSUGRA)-like unification of the soft parameters at the SUSY-breaking energy scale, such that 
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M = Mq 3 = M Ua = M Da = M L:i = M Ez = Msusy, 
M 1/2 = 2Mi = M 2 = ±M 3 , 
A = A t = A b = A T . 

These parameters are then fixed to their optimal values based on earlier studies [13} [T5] in order 
to minimise the set of scanned parameters. We then focus only on the effects of the Higgs sector 
parameters, which include the dimensionless Higgs boson couplings A and k along with their phases 
<j)\ and 4> K , as well as the soft SUSY-breaking parameters A\ and A K . From outside the Higgs sector, 
we only analyse the effect of the variation of the unified CPV phase of the third generation trilinear 
couplings, 4> Ao (= 4>A t = 4>A b =4>A T )- 

Before we discuss the three scenarios mentioned above, we note that the two heaviest Higgs 
boson mass eigenstates and always correspond to the interaction eigenstates and a in 
eq. Hence, a scenario is defined by the Higgs state that conforms to the LHC observations, out 
of the three light mass eigenstates, Hi, H2 and H3, and by the correspondence between the latter 
and the interaction eigenstates <p d , <p s and as- However, note that such a definition is adopted only 
so that a distinction between different scenarios can be made conveniently. Evidently, the behaviour 
of the 'observed' Higgs boson, H s [ g , with the CPV phases in a given scenario is a combined result 
of the set of parameters yielding that scenario rather than of its position among the mass-ordered 
Higgs states. The criteria for choosing the ranges of the scanned model parameters as well as the 
values of the non-Higgs-sector SUSY parameters thus depend on the scenario under consideration 
and are explained in the following. 

Scenario 1: In this scenario the lightest Higgs state, H\, is the SM-like one and corresponds to 4>% 
while H2 and H3 correspond to <p s and as, respectively. The requirement of obtaining a down- type 
Higgs state with mass close to 125 GeV and with SM-like couplings necessitates large soft SUSY 
masses and Aq. The values of /j, e s (= As, where s is the vev of S) and the gaugino masses are found 
to be in best agreement with the relic density constraints |15j . giving a neutralino with a large 
Higgsino component as the lightest SUSY particle. Further, A and k are chosen such that there is 
enough mixing of the doublet with the singlet Higgs boson so as to allow an H\ with the correct 
mass while keeping its couplings close to their SM values. We test two cases for this scenario, 
corresponding to two representative values of the parameter tan/3 (= v u /vd, where v u and Vd are 
the vevs of H u and H d , respectively), which is fixed to 6 in Case 1 and to 15 in Case 2. 

Scenario 2: This scenario is defined by the SM-like ~125 GeV Higgs boson being the second light- 
est Higgs boson, H2, of the model. There are two possibilities entailing such a scenario. It can 
correspond to 4> s in which case it has SM-like or bigger, as shown in |13j . We refer to this 
possibility as Case 1 of this scenario. It requires relatively large values of A and k, small values of 
the parameters A\ and A K and moderate values of soft SUSY-breaking parameters. H\ and H3 
in this case are (fi d and as dominated, respectively. For Case 2 of this scenario, we take a slightly 
different region of the parameter space which yields a H2 that is c^-like and hence has R~h 2 around 
the SM expectation. Hence, heavy unified soft squark mass and/or trilinear coupling are required 
in this case, but a light soft gaugino mass is preferred. A can be small to intermediate while k is 
always small. Finally, H\ and #3 are <f> s - and a^-like, respectively. 

Scenario 3: There also exists the possibility that the observed ~125 GeV Higgs boson is the H3 of 
the model which corresponds to (j) d , while both cj) s - and ag-like Higgs bosons are lighter. Such a 
scenario can be realised for very fine-tuned ranges of the parameters ^ an d A K for a given tan /3 

implying that after diagonalisation of the Higgs mass matrix, e.g., </)„ makes the dominant component of Hi. 
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Scenario 


1, Case 1 


1, Case 2 


2, Case 1 


2, Case 2 


3 


Fixed parameters 


M (TeV) 


5 


0.8 


3 


3 


M 1/2 (TeV) 


3 


0.35 


0.35 


1.5 


-A (TeV) 


10 


1 


4 


4 


ties (TeV) 


: 


L 


0.14 


0.14 


0.14 


tan /3 


6 


15 


1.9 


20 


10 


Scanned parameters 


A 


0.01 


- 0.1 


0.5 - 0.6 


0.01 - 0.3 


0.1 - 0.3 


K 


0.1 - 


- 0.3 


0.3 - 0.4 


0.01 - 0.1 


0.05 - 0.1 


A x (TeV) 


1.5 


- 3 


0.14 - 0.2 


0.2 - 0.6 


0.95 - 1.05 


-A, (TeV) 


1 - 


- 4 


0.2 - 0.25 


0.1 - 0.3 


0.07 - 0.09 



Table 1: Input parameters of the CPV NMSSM and their numerical values adopted in our analysis. 



value, with large soft squark and gaugino mass parameters preferred. Note that in this case the 
ag-like H% of Case 2 of Scenario 2 turns into H2 by becoming lighter than the c/>^-like state which, 
consequently, turns into H3. These two cases thus overlap slightly in terms of the relevant param- 
eter space of the model. 

We point out here that we do not consider a scenario with the ~125GeV Higgs boson corre- 
sponding to the as interaction eigenstate, since the pure (or nearly pure) pseudoscalar hypothesis 
is disfavored by the CMS Higgs boson analyses [36l[37j. Values of the fixed parameters as well as 
ranges of the variable parameters for all the above scenarios are given in Tab.HJ 

4 Scans and results 

We perform scans for each of the scenarios described earlier requiring the mass of -ff s i g (i.e., of H\ in 
Scenario 1, H2 in Scenario 2 and H% in Scenario 3) to lie in the range 124 GeV < rriH ni < 127 GeVH 
We additionally impose the condition Rff. > 0.5 on the signal Higgs boson. Furthermore, to each 
case in a given scenario corresponds four scans, such that in each of the scans only one of the 
following (set of) phases is varied: 

(i) 4> K , (ii) <j)x, (in) 4>a , (iv) 4>x and <p K , 

while fixing the rest of the CPV phases to 0°. Each scan thus checks the effect of a different CPV 
source at the tree-level and/or beyond. Note that, since the phases 4>x and 4> K enter at the tree-level, 
the contributions from the terms dependent on them reverse sign for (f>x, (ftn > 90°. This can lead 
to negative diagonal entries in the Higgs boson mass matrix, giving tachyonic mass eigenstates. 
On the other hand, 4>a does not have such a restriction and can yield physical masses even with 
values above 90°. However, for the sake of consistency, we only scan between 0° — 90° for all these 
phases. The relevant phase(s) is (are) varied in steps of 1° for scans (i)-(iii) above and in steps of 
5° for scan (iv). 

We should point out here that the measurements of the Electric Dipole Moment (EDM) of the 
electron, neutron and various atoms [38j EH HO] put constraints on the allowed values of 4>x and 

We thus use the central mass measurement of 125.5 GeV in accordance with the CMS results. 
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<f>A - However, the trilinear couplings of squarks and sleptons contribute to the EDMs only at the 
two-loop level and their phases are thus rather weakly constrained. One can, furthermore, assign 
very heavy soft masses to the sfermions of first two generations in order to minimize the effect of 
4>a on the EDMs, as pointed out in earlier studies for the MSSM [41]. In fact, such constraints 
can be neglected altogether by arguing that the phase combinations occurring in the EDMs can be 
different from the ones inducing Higgs boson mixing [42]. The phase of k, in contrast, has been 
found to be virtually unconstrained by the EDM measurements [32j[29]. 

Below we present our results separately for each of the five cases investigated. For evaluating 
the effect of the phases on w-# sig and R^j . qualitatively, we choose a set of four Representative 
Points (RPs), referred to as RP1, RP2, RP3 and RP4 in the following, for every scenario. RP1 
corresponds to a point for which the effect of at least one out of the three phases 4> K , and (f>A 
on rn,H sis is maximised for the given case. Similarly, RP2 is chosen such that the variation in R~ff. 
is maximal with one of the phases. RP3 and RP4 are points with the largest affect on rnH slg ana! 
R~ff. , respectively, observed in our scans with simultaneous variation of 4>\ and <p K . Note that in 
the discussion below the description of the behaviour of a given RP may not be equally applicable 
to all other good points, since it is chosen so as to understand the maximum possible impact of the 
phase that is the most influential one for a given case. In particular, for different RP2 and RP4 in 
a given case, the roles of 4>\ and 4> K could well be interchanged, since depends dominantly on 
the as component of Hi which is driven by a combination of A and k. 

4.1 Scenario 1: 

Case 1: In Fig. [TJ (top- left) we show the variation in the number of good points, i.e., points surviving 
the conditions imposed on rriH sig and Rj^ '. , for this scenario with small tan /3 for varying <j)\ , cj) K and 
4>a - Fig. Q] (top-right) shows the corresponding variation with (f>\ and 4> K changing simultaneously. 
The number of surviving points falls sharply with <j>A but relatively slowly with 4>\ and <p K . The 
variation is, however, not continuous since there are other parameters, A, k, A\ and A K , which are 
also scanned over for every value of a given phase. Moreover, the number of good points evidently 
depends on the conditions on rnH sig and R~ff. , so while both of these may be satisfied for one value 
of a phase one of these may be violated for the next. This is especially true for points producing 
7Bf/ sig and/or Rj^, near the boundaries of their required ranges. Note that the CPC case is subject 
to these two conditions also and, on account of being defined relative to this case, the number of 
good points does not represent all possible solutions for all values of the phases. Thus, it is likely 
that the CPC case for a given parameter set falls outside the defined ranges of "mH Ag and/or Rjj . , 
but the conditions on these are satisfied for a different value of a particular phase. Such a value 
of the phase can thus result in a considerable number of good points which were absent in the 
CPC case. Nevertheless, the aim here is to give an estimate of the effect of the CPV phases on the 
number of good points relative to the CPC case, rather than presenting a truly holistic picture. 
The values of other parameters corresponding to each RP for this Case are given below. 



Point 


A 


K 


A x (GeV) 


A K (GeV) 


RP1 


0.1 


0.11 


3000 


-1000 


RP2 


0.1 


0.25 


1500 


-3000 


RP3 


0.05 


0.11 


2000 


-4000 


RP4 


0.1 


0.25 


1500 


-3000 



Fig. Q] (middle-left) shows that the condition on m# si is the one mainly responsible for the drop 
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in the number of good points. In the figure we see m,H x falling slowly with increasing cj> K and <j>\, 
more so for the former than the latter. The drop in with increasing 4>a is relatively sharp. 
This is due to the fact that Hi in this scenario is very 2-like and the contribution from the singlet- 
like states is small since fairly small values of A and n are involved. Hence the phases of these 
parameters also have a lighter impact on itih 1 • On the other hand, in order to reach values up to 
125 GeV, the mass of Hi strongly relies on the trilinear coupling Aq and is consequently also more 
sensitive to its phase. Fig.Q] (middle-right) shows the combined effect of <p\ and <p K . We should 
point out again that at the tree-level only the difference 4>\ — (fi K enters, so a given combination of 
these phases can be mapped onto a particular value of cj)\ at the tree-level and should thus have 
the same effect. However, such a mapping is not possible when one includes higher order effects 
since at this level cp\ and (j) K appear independently of each other. The figure thus implies that the 
role of these phases is, expectedly, even smaller at the one-loop level. For small <j)\ larger values of 
4> K cause a bigger drop in itlh 1 while for larger 4>\ a smaller <p K results in more deviation from the 
CPC case. 

Bottom-left and -right panels in Fig. [I] correspond to RP2 and RP4, respectively, and show a 
fairly generic dependence of i?77. on the CPV phases for this scenario. We note that RTr. falls 
sharply with 4>\ and <j) K , more so for the former than the latter, until it reaches a minimum value 
and then rises again. The observed behaviour is due to the variation in the Hi — > bb width which 
increases sharply owing to the contribution from the as component of H\, for which the bb cou- 
pling gets enhanced faster than the scalar component. BR(iii — > bb) thus reaches a peak and then 
starts dropping. Expectedly Rj[. has negligible dependence on c/)a since the as component does 
not receive higher order corrections and thus does not have any dependence on Aq and its phase, 
while for the cpd component the variation in all BRs is almost constant when this phase is varied. 
RP4 in the figure on the right follows a trend similar to RP2, with R]j. showing greater variation 
compared to the CPC case for not too large values of and <p K . 

Case 2: For large tan (3 in this scenario the drop in the number of good points is even slower than 
in Case 1 as seen in Fig. [2] (top-left). For almost the entire range of <j) K all the points satisfy the 
conditions imposed, while in case of <p\ the number of good points falls quite abruptly around 75° 
and likewise for <f>A around 40°. A similar behaviour is observed in Fig. [2] (top-right) where the 
number of good points falls to for 4>\ > 75° irrespective of the value of cj) K . The RPs for this Case 
are as follows. 



Point 


A 


K 


A x (GeV) 


A K (GeV) 


RP1 


0.07 


0.11 


3000 


-1000 


RP2 


0.1 


0.25 


1500 


-3000 


RP3 


0.07 


0.11 


3000 


-1000 


RP4 


0.07 


0.28 


1833 


-3000 



The sudden fall in the number of good points is, once again, largely driven by TriH si „ which 
reaches its allowed lower limit gradually with increasing values of the phases as verified by Fig. [2] 
(middle-left and middle-right) for RP1 and RP3, respectively. The cut-off in the <p K line for RP1 
is due to the sudden drop of R]^. below 0.5. Such a behaviour of R]^ . is evident from Fig.[2] 
(bottom-left and bottom-right) for RP2 and RP4, respectively, where one sees that it remains al- 
most constant before dipping sharply for some given particular values of and <p K . This is again 
due to a sudden surge in the Hibb coupling and consequently in RR(Hi — > bb) which in turn causes 
BR(£/i —7- 77) to diminish, even though the Hi — > 77 width itself remains almost constant. Notice 
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Figure 2: Distributions of good points, rau aig and Rjj sig for Scenario 1, Case 2. 
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Scenario 


1, Case 1 


1, Case 2 


Points scanned for each <\> K or (p\ or 4>a q 


40000 


Points surviving for (j) K = = 4>a = 
mh Big for RP1 with 4> K = 4> x = 0a o = 
Kg. for RP2 with K = A = <p Ao = 


12494 
124.4570 
0.9990 


39976 
126.3754 
0.9997 


Max. points surviving, with <f) K 
Max/min m/j si obtained for RP1, with <j) R 
Max/min R 1 ^ . obtained for RP2, with (f> K 


12494, 
124.0058, 51 
0.9399, 10 


40000, 30-85 
124.8054, 88 
0.5001,64 


Max. points surviving, with <p\ 
Max/min . obtained for RP1, with cj)\ 
Max/ mm R H . obtained for KP2, with 0a 


12494, 
124.0009, 48 
0.7613,9 


39976, 0-4 
126.2300, 61 
0.9700, 68 


Max. points surviving, with 4>a 
Max/min m./j sig obtained for RP1, with <pA 
Max/mm it^ obtained for tirz, with c/>^ 


12494, 
124.0117, 19 
0.9990, 9 


39976, 0-3 
124.0015,41 
0.9998, 39 


Points scanned for each pair of <p K and 


10000 


Points surviving for (f) K = 4>\ = 4>a = 
m hs . g for RP3 with cj) K = (f) X = (pA = 
BCg s . for RP4 with K = A = Ao = 


3248 
124.4569 
0.9990 


9983 
126.3753 
0.9998 


Max. points surviving, with ^ k ,0a 
Min m^^ obtained for RP3, with <p K , <p\ 
Min R]j. obtained for RP4, with (p K ,4>\ 


3248, 0, 
124.0080, 75, 55 
0.7754,60,20 


10000, 5-10, 30-85 
125.797,5,85 
0.778, 115,40 



Table 2: Scan results for Scenario 1, Cases 1 and 2. All angles are in degrees. 



that since tan /3 is large in this Case, the hP d component of H\ has an enhanced coupling to bb while 
its as component is strongly dependent on A and k as well as their phases. The combined effect of 
these two factors is that sensitivity of the H\ — > bb decay to certain values of 4>\ and <p K is further 
enhanced through the relevant entry of the Higgs boson mixing matrix O. 

We should mention here that H2 and H3 are always very heavy, ~2TeV, for both the Cases of 
this scenario. Moreover, R?j Z follows the same trend with the variation in CPV phases as Rl7 
and is in fact always almost equal to it. Some important numbers corresponding to these Cases 
and to the four RPs in each of these are provided in Tab.[2j 

4.2 Scenario 2: 

Case 1: In Fig.[3] (top-left) we show the number of good points vs. the CPV phases. We see that 
the number of surviving points in the CPC case is much smaller compared to the two Cases of 
Scenario 1. However, in contrast to the latter, it falls more sharply with <fi\ and <p K than with 4>a - 
In fact, no good points are obtained for <ft\ > 25° when it is varied simultaneously with cfi K , as seen 
in Fig. [3] (top-right). This significantly different dependence of the signal Higgs boson in this case 
is due to the fact that it is now 05-like. In order to obtain such a singlet-dominated Higgs boson 
with the correct mass and 77 rate, fairly large values of the parameters A and k are needed, which 
result in an increased sensitivity to their phases also. On the other hand, since both (j) s and as 
states do not couple to matter, the only dependence on Aq and its phase comes through the small 
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doublet component of H 2 . The particulars of the RPs for this Case are given below. 



Point 


A 


K 


A x (GeV) 


A K (GeV) 


RP1 


0.51 


0.34 


187 


-239 


RP2 


0.52 


0.33 


200 


-233 


RP3 


0.578 


0.33 


160 


-233 


RP4 


0.6 


0.34 


187 


-211 



Fig. [3] (middle-left) for RP1 further verifies our statement above. In the figure we see that mH sis 
varies rapidly with (f)\ and 4> K but very slowly with 4>a - In fact, mH sis always lies within the allowed 
range for all values of <pA for RP1 which, however, is also due to the fact that fnu- is very close to 
the allowed upper limit in the CPC case. Notice also that fnjj^ first falls sharply with increasing 
4> K an d then rises again after reaching a minimum implying, again, that due to the large values of 
k and A, the composition of H 2 alters abruptly with these phases. H\ in this Case is </>[]-like and 
is mostly around 100 GeV while H% is as-like and around 300 GeV. Fig. [3] (middle-right) for RP3 
shows that the combined effect of simultaneously increasing <p\ and 4> K is that rnH sig increases very 
sharply with them. In fact, as in Scenario 1, the mixing effects cause the mass of Hi to fall with 
increasing values of the phases, which in turn results in H 2 gaining mass. 

Another characteristic feature of Case 1 of this Scenario is the fact that Rj[. can be fairly 
higher than the SM expectation and its dependence on 4>\ and 4> K is opposite to that on 4>a ■ As 
shown in Fig. [3] (bottom-left) for RP2, while i?77. increases with 4>a , it falls abruptly with (f)\ 

Slg 

and 4> K . The reason is that, while BR(H 2 — > 77) grows with both <j) K and 4>a , BR(H 2 — > gg), 
which is much more dominant, grows faster with <fi K than with (f)A - Further, BR(H 2 — > 66) remains 
almost constant with increasing values of the phases. In case of increasing <ft\, on the other hand, 
BR(H 2 — > 77) itself falls while both BR(H 2 — > gg) and BR(H 2 ->■ 66) grow. Evidently, with both 
4>\ and 4>k increasing simultaneously, Rl7. falls quite rapidly, as seen in Fig. [3] (bottom-right) for 

Slg 

RP4. Finally, R^ z in this Case is always slightly lower than FCQ ' . (e.g., it is ~0.9 for RP2) but 
shows a similar behaviour with varying CPV phases. 

Case 2: This Case corresponds to a i^-like H 2 state but shows a somewhat different behaviour from 
Scenario 1. In this Case the fastest but non-smooth variation in the number of good points comes 
from <j) K , as seen in the top- left panel of Fig.[U Also, in contrast to Scenario 1, the number of good 
points, which is fairly low for the CPC case, rises slowly with increasing <pA and then starts falling. 
In the top-right panel of the figure we see that when both <fi\ and <p K are varied simultaneously, the 
drop in the number of good points is very slow, much like it is when (f>\ is varied alone (left panel). 
Below we give the values of other Higgs sector parameters for the four RPs of this Case. 



Point 


A 


K 


A x (GeV) 


A K (GeV) 


RP1 


0.047 


0.14 


600 


-180 


RP2 


0.044 


0.015 


467 


-147 


RP3 


0.04 


0.012 


556 


-153 


RP4 


0.05 


0.014 


289 


-160 



The middle- left panel of Fig. U for RP1 shows that ™-H sig rises continuously with increasing 4> K 
and but falls with (pA - The increase in mass of H 2 due to 4>\ and 4> K is, in analogy with 
Case 1, due to the mixing effects with the as-like H%, whose mass falls with these phases. On the 
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other hand, since now H2 itself is </>^-like, the dependence on Aq and consequently on its phase is 
enhanced again, compared to the Case 1. We note here that m# 3 for the CPC case of RP1 is fairly 
small, ~ 180 GeV, and falls abruptly with <j) K but rises slowly with <p\. The middle-right panel for 
RP3 shows that the simultaneous variation in <f)\ and (j) K also results in a slow increase in the mass 
of H2. However, the difference in mass is maximised compared to the CPC case when one of the 
two phases has a small value while the other has a large value. 

In the bottom panels of Fig. H] we show the variation in R~Q with the CPV phases for RP2 

and RP4 of this case. When CP is conserved in this Case . is generally slightly lower than 1 
(SM expectation) but sits well within the observed range reported in the latest results by the CMS 
collaboration [5]. Moreover, the variation with the CPV phases is comparatively slower than the 
Cases discussed so far, with the strongest dependence b eing on (f) K . Rlfj . falls smoothly with (pAo 
but with increasing cj)\ it first rises slightly and then starts falling slowly. This is due to the fact 
that BR(if2 — > 77) itself falls with increasing (J)a and <p K but rises slowly with In addition, 
BR(#2 - > gg) also falls with and 4> K while BR(i?2 — > bb) remains almost constant for these 
phases. However, BR(#2 — > ZZ) rises with <f)\ causing RQ . to fall beyond a certain value of the 

sig 

latter. In case of 4>a q BR(i?2 — > ZZ/WW) falls rapidly while BR(.£f2 — > bb), which is the most 
dominant of all BRs, rises faster, causing in turn iCff. to drop. For RP4 in the bottom-right panel, 
the combined effect of <j>\ and (f) K is that ICQ. falls slowly, except for small values of cj>\ and large 

sig 

values of <p K when it tends to rise but only minimally so. Once again, Rfj Z is very close to ICQ. 
for RP2 and RP4 in this Case also and follows a similar trend in variation with an increase in any 
of the three CPV phases. 

Some particular values corresponding to the benchmark points for the two Cases of this Scenario 
are given in Tab.[3l 

4.3 Scenario 3: 

Although, as noted earlier, the parameter space corresponding to this Scenario overlaps a little 
with the Case 2 of Scenario 2, some significant differences are noticeable. According to the top-left 
panel of Fig. [5] a considerably larger number of points survives the scans for the CPC case. The 
number of good points falls slowly with cpA and , more so for the former than the latter, but 
with increasing <j) K it first rises slightly and then starts falling very slowly. When (f>\ and (j) K are 
increased simultaneously, the drop is sharper than with the two phases varied individually, as seen 
in the top-right panel. The representative points of this Scenario have the following co-ordinates. 



Point 


A 


K 


A x (GeV) 


A K (GeV) 


RP1 


0.216 


0.09 


983 


-76.7 


RP2 


0.163 


0.055 


950 


-70 


RP3 


0.189 


0.05 


983 


-85.6 


RP4 


0.167 


0.056 


983 


-87.8 



The middle- left panel of Fig. [5] shows that, similarly to the Case 2 of Scenario 2, mH sig increase 
with and (j) K (although now the former has the largest effect of the two) while it falls with 
increasing (pA - However, as noted earlier, such a behaviour due to <j)\ and <p K is specific to this 
particular RP1 and may not be the same for a different RP1. In the middle-right panel of Fig. [5] 
for RP3 we see that, as with the corresponding point of Scenario 2, mH aiK rises as and (f> K 
increase simultaneously, showing maximum enhancement when both the phases have large values. 
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Scenario 


2, Case 1 


2, Case 2 


3 


Points scanned for each cj> K or 4>\ or 4>a 


40000 


Points surviving for 4> K = 4>\ = (j)A = 
m hs . g for RP1 with (p K = cpx = 4>A = 
R^ s . s for RP2 with c/> K = c/> A = cj) Ao = 


6017 
126.8907 
1.25 


8792 
125.1155 
0.8551 


16701 
125.9961 
0.8701 


Max. points surviving, with K 
Max/min m,h sig obtained for RP1, with (p K 
Max/min Rj[. obtained for RP2, with 4> K 


6017,0 
125.7728, 18 
1.2413, 13 


8792, 
126.7120,59 
0.8430, 16 


17011,22 
126.9878, 31 
0.8294, 60 


Max. points surviving, with 

Max/min m^. obtained for RP1, with <p\ 

Max/min R]^, obtained for RP2, with cf>\ 


6017,0 
124.1719, 11 
1.2253,8 


8792, 
125.5718,59 
0.8419, 57 


16701,0 
126.9616, 16 
0.8009, 39 


Max. points surviving, with <j)A 
Max/min m^. obtained for RP1, with 4>A 
Max/mm Rjj . obtained for RP2, with <pa 


6017,0 
124.2631,88 
1.3067, 42 


9126, 16 
124.0073, 28 
0.8280, 45 


16701,0 
124.0672, 44 
0.8683, 17 


Points scanned for each pair of (f) K and <j)\ 


10000 


Points surviving for (fi K = = 4>a = 
m hs . g for RP3 with (p K = cpx = 4>A = 
R% s . r for RP4 with (j) K = (j) X = cj) Ao = 


1555 
124.3536 
1.0609 


2141 
124.8334 
0.8851 


3862 
124.2729 
0.8899 


Max. points surviving, with (j) K ,4>\ 
Max m/i sig obtained for RP3, with (p K , 
Min Ry , obtained for RP4, with <ft K ,(J)\ 

sis 


1555,0,0 
126.8948, 25, 15 
0.5385,20,30 


2141,0,0 
126.4176, 55,55 
0.8426,55,30 


3862,0,0 
126.9995,40,45 
0.5015,20,75 



Table 3: Scan results for Scenario 2, Cases 1 and 2, and Scenario 3. All angles are in degrees. 



We should point out here that in this Case the H\ (almost always i^g-like) is always lower than 
~115GeV and can be as low as 20GeV while H2 (mostly as-like) ranges between 70GeV and 
124 GeV (the upper limit is explicitly imposed in the scans). 

In the bottom-left panel of Fig.0 R^ff. for RP2 is shown as a function of the various phases. 
Once again, in analogy with Case 2 of Scenario 2, the value of Rjj. is slightly lower than the SM 
expectation. While it remains almost constant with 4>a , it falls slightly more sharply with (j>\ 
than with (j) K . In case of increasing <pA BR(i?3 — > 77) itself falls slowly while BR(//3 —> bb) rises. 
With increasing cf)\, on the other hand, while BR(i^3 — > gg) and BR(i^3 — > bb) both fall slowly, 
BR(i^3 — > WW/ZZ) rises comparatively rapidly. A similar pattern is seen for <f> K also, although 
there BR(i^3 — > 77) itself rises minimally with increasing phase, hence resulting in a slower drop 
in Rj[ . compared to the case for <p\. In the bottom-right panel of Fig. [5] we show the effect on 

sig 

R~H . of simultaneously changing and <p K for RP4 of this Scenario. We note that for large <j)\ 

sig 

and small 4> K Rjj '. can be considerably lower than its value corresponding to the CPC case. 

sig 

As with all the other Cases with a c^-like H sig , Rfj Z has similar values as R]^. when the CPV 
phases are zero and a similar behaviour when these phases are varied. Finally, some details relevant 
to the four RPs of this Scenario are given in Tab.O 
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Figure 3: Distributions of good points, niH sig and Rjj . for Scenario 2, Casel. 
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Figure 4: Distributions of good points, Tau aig and Rjj sig for Scenario 2, Case 2. 
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5 Summary 



In summary, we have demonstrated that the CPV NMSSM offers some interesting solutions to 
the LHC Higgs boson data, which differ substantially from well-known configurations of the CPC 
NMSSM, thereby augmenting the regions of parameter space which can be scrutinised at the CERN 
collider. We have concentrated on the case in which only three CPV phases, 4> K , 4>\ and 4>a , enter 
the Higgs sector. We have then checked the two-fold impact of these phases, varied separately or 
simultaneously but always independently from each other, on the mass as well as signal strength 
of the assumed signal Higgs boson in the 77 decay mode, in different model configurations. 

The overall picture that emerges is that any of the three lightest Higgs states of the CPV 
NMSSM can be the one discovered at the LHC. We have illustrated this by using five benchmark 
cases in the parameter space of the model that can easily be adopted for experimental analyses. Our 
analysis also proves that the possibility of explicitly invoking CPV-phases is not ruled out by the 
current LHC Higgs boson data in any of our tested plausible NMSSM scenarios. Finally, a numerical 
tool for analysing the Higgs sector of the CPV NMSSM has also been produced and is available 
upon request. Obvious outlook of this analysis will be to consider the possibility that companion 
Higgs boson signals to the one extracted at the LHC may emerge in the CPV NMSSM, so as to 
put the LHC collaborations in the position of confirming or disproving this SUSY hypothesis. An 
investigation on these lines is now in progress. 
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A Mass Matrices 

Detailed expressions for the one- loop Higgs boson mass matrices can be found in [28j 123 ED]- 
Here we only reproduce the tree-level mass matrix to show the dependence on <fi\ and <p K since 
the dominant contributions from these phases arise at this level. Note that the tree-level sfermion, 
neutralino and chargino mass matrices given below are complex by definition. The one-loop effective 
Higgs potential receives further contributions from (f>\ and (f>A Q through the squark and stau sectors 
and from and <j) K through the neutralino and chargino sectors. 

• The neutral Higgs boson mass matrix may be written as: 



Using the minimisation conditions of the Higgs potential, one can define some convenient 




(A.l) 



parameters: 




(A.2) 
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with 

<j)' x = <t>\ + + y and <// K = <j) K + 3<p . (A.3) 

In terms of these parameters, the entries of the top left 3x3 CP-even block in eq. (lA.lh are 
given as 





92 + ST 
4 




92 + gl 
4 




D VdVu 

= Kx 

vs 


(M 2 S ) 12 


= (MD21 


(Ml)l3 


= (M 2 S ) 31 


(-MD23 


= (M 2 S ) 32 



rf,+ [R x + ±Kv s 



< + [R\ + z nv s 

+ 2\k\ 2 v 2 s - R K v s , 



v u vs 

Vd 

Vdvs 



gl + g\ 



1 



+ \M )v d v u - (Rx + -nvs)v s , 



(A.4) 



The bottom right 2x2 CP-odd block reads 



M 2 P 



I'S 



(R\ - lZv s )v 



(R x -Kvs)v Rx^+2Kv d v u + 3R K v s 



(A.5) 



Finally, the entries of the off-diagonal block responsible for mixing between CP-even and 
CP-odd states is given as 

\Lv u v s 

-\lvdV S I • (A.6) 



M 2 SP 








-\Xv u vs -^Tvdvs 



The chargino mass matrix, in the (W ,H ) basis, using the convention H L , R s = H d , u y can 
be written as 



M 



c 



M 2 V2M W cos /3 



(A.7) 



which is diagonalised by two different unitary matrices as CrM.qC\ = diag{m~±, m~±}, 
where m~± < m~±. 

Al A 2 

The neutralino mass matrix, in the (B, W , H®, S) basis, can be written as 

/Mi -M z cos Ps w M z smj3s w \ 

M 2 M z cos j3c w -M z sin /3c w 






M 



N 



\ 



V2\k\v s J 



(A. 



It is diagonalised as N*MnN^ = diag {rn~o , m^o , m~o , m~o , m^o ) , where N is a unitary matrix 



and m^o < m^o < m^o < rrir-o < m~o. 

Xl — A2 — A3 — A4 — A5 
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For the stop, sbottom and stau matrices, in the (qL,QR) basis, we have 

~ 2 _ / M| + m 2 + cos2/3Mf (§-§ S ^) h* t v u (A* u e~ M - ^e^+^ cot (3)/V2 \ 
M ' ~ 1 ^(Xe^-^e-^+^cot^/v^ M| + ro? + cos 2 (3 M 2 Q t s 2 w J ' 



M| 3 + m 2 + cos2/?Mf (-§ + \s 2 w ) h* b v d (A* b - J^ e *(^+^) tan/3)/v^ \ 
- J^e- i (^ +0 +^) tan 0)/V2 M| g + + cos 2/?Mf Q b s 2 w J 



Mi 



+m 2 T + cos2/3M| (s 2 w - 1/2) - » e ^ +e+ ^ tan/3)/\/2 \ 

E 3 



h T v d (A T - l^. e -<4>x+<>+<p) tan p) /y/2 + m 2 T - cos 2/3M§ s 2 w ) 



where yt, yb and y T are the Yukawa couplings of t, 6 quarks and r lepton, respectively, and 
s\y = sin^vF- These sfermion mass matrices are diagonalised by unitary matrices U* such 
that TJft M. 2 f Uf = diag(m~ , m~ ) with m~ < m~ for / = t, b and r. 

7 /l /2 /l /2 
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